Spectral Degeneracies in the Totally 
Asymmetric Exclusion Process 



O. Golinelli, K. Mallick 



Service de Physique Theorique, Cea Saclay, 91191 Gif, France 

December 16, 2004 



Preprint T04/167; .arXiv:cond-mat/0412462| 



Abstract 

We study the spectrum of the Markov matrix of the totally 
asymmetric exclusion process (TASEP) on a one-dimensional pe- 
riodic lattice at arbitrary filling. Although the system does not 
possess obvious symmetries except translation invariance, the 
spectrum presents many multiplets with degeneracies of high or- 
der. This behaviour is explained by a hidden symmetry property 
of the Bethe Ansatz. Combinatorial formulae for the orders of 
degeneracy and the corresponding number of multiplets are de- 
rived and compared with numerical results obtained from exact 
diagonalisation of small size systems. This unexpected structure 
of the TASEP spectrum suggests the existence of an underlying 
large invariance group. 
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1 Introduction 

The asymmetric simple exclusion process (ASEP) is a driven lattice gas 
model in which particles interact by hard core exclusion. This simple system 
has been introduced as a building block for models of hopping conductivity, 
motion of RNA templates, traffic flow and surface growth. From a theo- 
retical point of view, the ASEP plays a fundamental role in the study of 
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non-equilibrium processes: many exact results have been obtained concern- 
ing one-dimensional phase transitions (Derrida et al. 1993), phase segrega- 
tion (Evans et al. 1998), large deviations functions and fluctuations far from 
equilibrium (Derrida et al. 2003). For a review, see Derrida (1998), Schiitz 
(2001). In the absence of a driving field, the symmetric exclusion process 
can be mapped into the Heisenberg spin chain. The asymmetry due to a 
non-zero external driving field breaks the left/right symmetry and the ASEP 
is equivalent to a non-Hermitian spin chain of the XXZ type. The ASEP 
can also be mapped into a two-dimensional six- vertex model at equilibrium. 
These mappings allow to use the methods of integrable systems, such as the 
Bethe Ansatz (Dhar 1987, Gwa and Spohn 1992, Kim 1995). 

In a recent article (Golinelli and MaUick 2004), we carried out a spectral 
study of the Markov matrix of the exclusion process on a periodic lattice. At 
half-filling the system is invariant under charge conjugation combined with 
rcficction in addition to being translation invariant. We showed that these 
symmetries predict the existence of singlets and doublets in the spectrum. 
However, for the totally asymmetric simple exclusion process (TASEP) the 
spectral structure is much richer. We observed numerically that unexpected 
degeneracies of higher order exist and become generic as the system size 
increases. We explained, in an heuristic manner, the existence of these de- 
generacies by using the fact that some of the solutions of the Bethe equations 
appear in pairs with opposite values. 

In the present work, we generalize our previous study to the TASEP at 
arbitrary filling. Wc perform an exhaustive analysis of the spectral degenera- 
cies of the TASEP on a periodic lattice. The orders of degeneracies observed 
and the number of multiplets with a given degeneracy depend on commen- 
sur ability relations between the number of sites and the number of particles. 
Although the Bethe equations do not exhibit any obvious symmetry, we find 
that they possess an invariance under exchange of roots that allows to predict 
combinatorial formulae for the orders of degeneracies and the number of mul- 
tiplets of a given order of degeneracy. Our formulae are confirmed by direct 
numerical diagonahsation of small size systems. This peculiar structure of 
the TASEP spectrum suggests the existence of some underlying symmetries 
of the model that may shed light on its remarkable combinatorial properties. 

The outline of this article is as follows. In section 2, the definition and 
basic properties of the TASEP are recalled. In section 3, we present a self- 
contained derivation of the Bethe equations based upon the fact that the 
Bethe wave function is a determinant for the TASEP model. The symmetries 
of the Bethe equations are studied in section 4 and combinatorial expressions 
for the degeneracies arc derived in section 5. In section 6, wc discuss the 
behaviour of large size systems and present numerical results. Concluding 
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remarks appear in section 7. In the Appendix, the geometrical setting of the 
roots in the complex plane is briefly described. 



The simple exclusion process is a stochastic process in which particles hop 
on a lattice and respect the exclusion rule that forbids two or more particles 
per site. On a one-dimensional lattice, this rule prohibits overtaking between 
particles. 

In this article we consider a periodic 1-d lattice of length L, i.e., a ring 

where sites i and i + L are identical. The system is closed and the number N 
of particles is fixed with N < L. The filling (or density) is given by p = N/ L. 

The particles evolve with the following dynamics rule: during the time 
interval [t,t + dt], a particle on a site i jumps with probability dt to the 
neighbouring site i + 1 if it is vacant. As the jumps are allowed in only one 
direction, the model considered is the totally asymmetric exclusion process 



A configuration C of the system is characterised by the list of the N occu- 
pied sites amongst the L available sites. The total number of configurations 
is therefore 



Let iptiC) be the probability that the configuration of the system at time t 
is C. As the TASEP is a continuous-time Markov (i.e., mcmoryless) process, 
the fi-dimensional vector tpt evolves according to the master equation 



where M is the ^7x^2 Markov matrix. For C ^ C, the term M{C',C) 
represents the transition rate from C to C: it is equal to 1 if C is obtained 
from C by an allowed jump of one particle, and is otherwise. The element 
—M{C, C) is equal to the number of allowed jumps from C. Thus, the sums 
columns of the Markov matrix vanish and the total probability is conserved. 
For the exclusion process on a periodic lattice, the sums over lines of M also 
vanish: the stationary probability (which corresponds to the eigenvalue 0) is 
thus uniform: ip{C) = 1/Vt. 

As the dynamics is ergodic and aperiodic, M satisfies the conditions of 
the Perron-Probenius theorem: the eigenvalue is non-degenerate and all 
the other eigenvalues have strictly negative real parts (equal to the inverse of 



2 The TASEP model 



(TASEP). 




(1) 
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the relaxation times). As M is a real non-symmetric matrix, the eigenvalues 
are either real numbers or complex conjugate pairs. 

In this work, we shall investigate the spectral degeneracies, i.e., equalities 
amongst the eigenvalues of the Markov matrix. 



3 Derivation of the Bethe equations 

Since the work of Dhar (1987), it is known that the Bethe Ansatz can be 
applied to the ASEP. In this section, we give a self-contained derivation of 
the Bethe equations for the particular case of the TASEP, much simpler than 
that of the generic ASEP (Gwa and Spohn 1992). 

A configuration C will be represented by the sequence (xi,X2, ■ ■ ■ ,xn), 
the integers Xj being the positions of the particles with 

1 < Xi < X2 < ■ . ■ < Xjy < L . (3) 

The idea of the Bethe Ansatz consists in writing the eigenvectors ip of the 
Markov matrix as linear combinations of plane waves (see, e.g., Gaudin 1983). 
In fact, the Bethe wave function ip of the TASEP is a determinant (Gaudin 
and Pasquier 2004). We therefore define ip as 

7/'(xi,...,X7v) = det(i?), (4) 

where R is a N x N matrix with elements 

R{z,j) = for 1 < < A^, (5) 

{zi, . . . , zn) being given complex numbers. If we assume that ip is of this 
form and that it is an eigenvector of M, the Zj's then must satisfy some 
conditions, the Bethe equations, that we now re-derive. 

We first show that ip defined by equations (jH Ej) satisfies two identities 
which are valid for any values of Zi and of Xj, even without imposing the 
ordering given in equation Q. The first identity is 

N 

Eip{xi, . . .,xn) = '^bP{xi, ...,Xk-l,.. .,xn) - ip{xi, ...,Xk,.. .,Xn)] , 

k=l 

(6) 

for any {xi, . . . ,xn) and {zi, . . . , zn), and where E is given by 

N 

E = -N + J2l/z,. (7) 

1=1 
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Equation (0) is obtained by writing 

1p{xi, . . . ,Xk - I, . . . , xn) - ^(Xi, ...,Xk,..., Xn) = 

det (^R{i, 1), . . . , (^1 - R(t, k),..., R{t, N) ) . (8) 

This determinant is similar to det(i?) except for the k-th column. Expanding 
this determinant over all permutations of {1, ... , A^} and performing the sum 
over k = 1, . . . , N leads to equations (0 EI)- 

The second identity valid for any {zi, . . . , zjsi) and any (xi, . . . , xjs) with 
Xk-i = Xk (two particles collision), is 

1p{Xi, ...,Xk,Xk,...,Xn)- ^{Xi, ...,Xk,Xk + l,...,Xn)=0. (9) 

The left hand side of equation (0) can be written as det(i?) where i? is a 
matrix that is identical to R but for its k-th column that is given by 

^) = '^iiZ'^^ = (1 _t])fc-i = ^(^' ^ - 1) = ^(^' ^ - 1) (10) 

The (A;— l)-th and the A;-th columns of R are equal and, therefore, det(i?) = 0. 
This proves equation Q. 

The eigenvalue equation, Eip = Mip, is written as equation ^ except 
that the sum is restricted to the allowed jumps of particles, i.e., to the indices 
k such that Xk-i + 1 < Xk- However, in equation ©, the terms with Xk-i + 1 = 
Xk vanish thanks to equation @. Thus equation © is identical to the 
eigenvalue equation if the eigenvector has the form assumed in equations (jH 

ED- 

The function ip must also satisfy periodic boundary conditions 

iIj{Xi,X2, ...,Xn)= 1P{X2, . . . , X^, Xi + L) . (11) 

The periodic conditions are the ones that quantify the eigenvalues by impos- 
ing a set of equations on the Zi's, the Bethe equations. Denoting by i and j 
the generic line and column indices of the matrices, respectively, we have 

i.ix2,...,x.,x^ + L)=det^-^,...,j^,...,^^^j . (12) 

By cyclic permutation of the columns, we obtain 

^(X2, . . . ,Xn,Xi + L) 
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If-' det 



= i-^f-'Ili^-^k) det(^^^-^i?(z,l),...,i?(z,j),...j .(13) 

This last term is equal to 'iIj{xi,X2, ■ ■ ■ , x„) = det(i?) if Zi, . . . , zn are solu- 
tions of the N Bethe equations 

N 

(^,-l)V^ = -n(l-^fc) for 1 = 1,. ..,N. (14) 

k=l 

The vector ip defined by equations (j3J EI) is then an eigenvector of the Markov 
matrix M with eigenvalue E given by equation (|7j). 

The uniform stationary probability with E = corresponds to the very 
special solution where all the Zi = 1. For all the other solutions of the Bethe 
equations, the z^'s are distinct and are different from 1; hence the determinant 
det(i?) does not vanish. 



4 Symmetries of the Bethe equations 

In this section we show that the Bethe equations p4|l have certain solutions 
that are distinct but lead to the same eigenvalue E. 

Following Gwa and Spohn (1992), we introduce new variables Zi = 2/zj — 1 
in the Bethe equations (|Tl|l which then become 

(1 - f.)^ (1 + f.)^-^ = -2^ n ^ for t = l,...,N. (15) 

k=l i- 

The corresponding eigenvalue E is given by 

N 

2E = -N + J2 4 • (16) 

k=l 

We remark that the right hand side of equation ()15|) is independent of i. We 
shall analyse equation (fTH|l in three steps. Firstly, we consider the polynomial 
equation of degree L for a given complex parameter Y, 

{1- Zf {1 + Z)^-^ = Y . (17) 

We call {Zi, . . . , Zl) the roots of this polynomial. In the Appendix, we 
explain how the ZiS can be labelled so that for a given i, Zi is an analytic 
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function of Y in the complex plane with a branch cut along the real semi-axis 
[0, +cxd). Secondly, we choose a set c = {ci, . . . ,C]y} of distinct indices 
among {1, . . . ,L}. The number of possible sets c is precisely Q, the total 
number of configurations (equation (Q)). Finally, for a given choice set c, we 
define a function of Y 

^7-1 

A,(F) = -2^ n • (18) 

k=l ' 

The Bethe equations (|TH|l are now equivalent to the self-consistency equation 

A,{Y)=Y. (19) 

For a given choice set c and a root Yc of the last equation, the Z^s are deter- 
mined by equation (fTTj) . The solutions of the Bethe equations are then given 
by Zk = ^Cfc- The corresponding eigenvalue is obtained from equation (fTBj) 

2E, = -N + Y.Z,^. (20) 

k=l 

The eigenvector ip is given by equations (jH El), using Zi = 2/{zi + l). 

In order to understand the origin of the spectral degeneracies, we must 
consider the case where L (the number of sites) and N (the number of par- 
ticles) are not relatively prime. We define the integers p, n and / as follows 

p = gcd(L,A^), L = pl, N = pn. (21) 

Let (?/i, . . . , Hp) be the p-th roots of Y (i.e., yl = Y) labelled as < arg(yi) < 
. . . < arg(?/p) < 271. Equation (fT7|) is thus equivalent to the p polynomial 
equations of degree / 

Qfc(Z) = (l-Z)"(l + Z)'-"-y, = for k = l,...,p. (22) 

Thus the set {Zi, . . . , Z^} of the L solutions of equation (fTTjl is made up of p 
packages, the k-th package being constituted by the I solutions of Qk{Z) = 0. 
Let us call Vk the set of indices of the fc-th package: in other words the 
solutions of Qk{Z) = are the Zi with i G Vk- The labelling of the Zi 
described in the Appendix shows explicitly that 

Vk = {k,k + p,k + 2p,...,k + L-p} , (23) 

i.e., Vk contains the indices i such that i = k modulo p. 
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For any given package Vk, we have the fundamental equations 



2n-l, 



(24) 



n 



Zi + l 



1. 



(25) 



We emphasize that the right hand sides are independent of k. These identities 
are derived as follows. We have {—l)"'Qk{Z) = Y[i(^Vki^ ~ ^i) because the 
Zj's with i EVk are the roots of the polynomial Qk{Z). The evaluation of the 
coefficient of Z^^^ leads to the equation except when / = 1. Similarly, 
the evaluation of QkiX) IQki,—^) yields equation (j25|l except when n = 
or n = /. These exceptions correspond to a trivial model which is either 
empty, iV = 0, or full, N = L, (the spectrum of M is then reduced to the 
single eigenvalue {0}). In the following, we assume that < N < L and so 
equations EH)) are true. 

We now consider a solution Yc of the Bethe equation (|T9|l associated with 
a given choice set c. Assume that there exists a package Vf such that c 
contains Vf (i.e., Vf C c) and that there exists another package Ve such that 
c is disjoint from Ve (i.e., H c = 0). We define a new choice set c' of N 
indices by exchanging the package Vf with Ve, i.e., c' = (c/Vf) UVe- We now 
show that the eigenvalues associated with c and c' are equal. Indeed, because 
of equation the contribution of Vf in AdY) and the contribution of Ve in 
Ac'(Y) are both equal to 1 and therefore Ac(Y) = Ac'(Y), the other packages 
contained in c and c' being the same. Thus Yc is also a solution of the Bethe 
equation associated with the set c', i.e., Ac'iYc) = Yc. Besides, thanks to 
equation (j^ . we notice that the contribution of Vf to the eigenvalue Ec is 
equal to the contribution of Ve to Ec'. Thus, because the other packages 
contained in c and c' are the same, we conclude from equation ^K)^ that 
Ec = Ec'. However, the corresponding eigenvectors are different: some of 
the chosen Zi^s being different for c and c', the functions ipd^i, . . . ,xn) and 
ipc'i^i, . . . , xat) are not equal for the fl different configurations. We have thus 
obtained a degenerate eigenvalue Ec = Ec' associated with two different sets 
c and c'. 

To summarise, an eigenvalue corresponding to a choice set c is degenerate 
if there exists at least one package Vf entirely contained in c and at least one 
package Ve that does not intersect with c. The fundamental reason is that a 
full package of Zj's does not contribute to the Bethe equations and adds up 
to a constant contribution in the eigenvalue formula. Therefore exchanging 
the packages Vf and Ve leads to the same eigenvalue but not to the same 
eigenvector and results in degeneracies in the spectrum. 
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As the size of the packages is I, we note that degeneracies are possible 
only if / < < L - / with / = L/ gcd(L, A^). 

5 Combinatorial formulae for the degenera- 
cies 

We shall now enumerate the degeneracies in the spectrum of the Markov 
matrix M of the TASEP with A^ particles evolving on a ring of L sites. We 
assume the following one-to-one hypothesis: for each choice set c (among the 
fl possible sets), the self-consistency equation flT^ has a unique solution 
that provides the eigenvalue Ec and the eigenvector tp^- We further assume 
that these eigenvectors are linearly independent; this hypothesis, combined 
with the fact that the dimension of the configuration space is Q, implies that 
the Bethe equations provide a complete basis of eigenvectors. We emphasize 
that the one-to-one hypothesis is stronger than the assumption that the Bethe 
basis is complete. We have observed numerically on small size systems that 
the functions Ac{Y) are usually contraction mappings (which would imply the 
one-to-one hypothesis), but we have not succeeded yet to obtain a rigorous 
proof. However we will see that this one-to-one hypothesis allows to count 
the degeneracies and that the results are in perfect agreement with numerical 
diagonalisations. Thus we are convinced that this hypothesis, or a weaker 
hypothesis with the same consequences, is true and that it should be possible 
to prove it. 

With this one-to-one hypothesis, counting degeneracies becomes merely 
an exercise in combinatorics. We first introduce some notations. We recall 
(see Section H} that c = {ci, . . . , c^} is a set of A^ integers chosen amongst 
{1, . . . , L}; moreover {1, . . . , L} is partitioned in p packages Pi, . . . , Vp, each 
containing I = L/p integers with p = gcd(L, A^) (see equat ions (fTTH^ ) . For 
a given set c and for < z < /, we denote by the number of packages Vk 
with i elements in c (i.e., such that Vk (~] c has i elements). Thus oq is the 
number of packages that do not intersect c: such packages will be referred to 
as 'empty' packages. The number of the 'full' packages (i.e., entirely included 
in c) is ai. We call partial packages those that are neither empty nor full. 
Following this definition, we have 

/ I 

ai>0 , ^ai= p , ^iai = N , (26) 

1=0 i=0 

p being the total number of packages and A^ the cardinality of the set c. We 
call a = (oq, ai, . . . ,ai) an admissible partition if it satisfies the relations in 
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equation Equivalently an admissible partition corresponds to a parti- 

tion of the integer in which each term is < / and which contains at most 
p terms. 

The total number uj{a) of choice sets c corresponding to a given admissible 
partition a is 



o) (l) (2) -U -^'S^U .(27) 



=0 



In this equation, the first factor enumerates the number of choices for each 
type of packages among the p available packages. The factors of the type 
[J give the number of choices of i elements among / for each of the Oj 
packages. According to the one-to-one hypothesis, Lj{a) represents also the 
number of eigenvalues associated with the admissible partition a. 

We have shown at the end of the previous section that two sets c and 
c' provide a degenerate eigenvalue Ec = Ec' if they are built from the same 
partial packages and differ only by the selected empty and full packages. 
Thus the eigenvalue E^ is d{a) times degenerate with 

d(a)=(^^° + ^'). (28) 

This relation is obtained by enumerating all the choice sets c' obtained from 
c by keeping the partial packages unchanged and choosing ai full packages 
from the remaining ai + packages. We remark that c and c' correspond 
to the same admissible partition a. We also emphasize that the degeneracy 
order depends only on a and not on the precise choice set c. To resume, a 
single order of degeneracy d{a) is associated with the admissible partition a. 

Consequently, the u^a) eigenvalues corresponding to the admissible par- 
tition a form m(a) multiplets of d{a) degenerate eigenvalues, where m(a) is 
given by 

d{a)- a,W...ai.,\{ao + a,y\l) [2) ■■■[l-lj "^^^^ 

Because the value of E^ depends only on the roots Zi belonging to the partial 
packages, we remark that this equation can also be obtained by enumerating 
the number of choices for these roots: the first factor counts the number 
of choices for the partial packages among the p packages and the other fac- 
tors enumerate the choices of i elements among / for each of the Oj partial 
packages. 
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In order to know the total number ■m{d) of multiplets with degeneracy of 
order d, we must sum over all admissible partitions a with d{a) = d, 

m{d) = J2 ^(«)- (30) 

a; d{a)=d 

In the particular case of half- filling {L = 2N), we have p = N and 1 = 2, 
n = 1. The admissible partitions and the corresponding degeneracies are 
parameterised by oq: equation (j26|l leads to a2 = ao, ai = N — 2ao and 
d{a) = (^^°^- The relation between the admissible partitions and orders of 
degeneracies is therefore one-to-one (i.e., two different admissible partitions 
have different orders of degeneracies): thus the sum in equation (j30|l reduces 
to a single term. Explicit formulae and numerical results are given in (Go- 
linelli and Mallick 2004). 

For fillings other than 1/2, we have I > 2 and two different admissible 
partitions can lead to the same order of degeneracy. Equation ()30|) can 
not be further simplified and the enumeration of admissible partitions seems 
mandatory. Nevertheless, we can verify the sum rule 

n = Y.d m{d) = Y: d{a) m{a) = E^'' fl A i^T (31) 

d>l a a i=0 V/ 

where the last sum runs over the admissible partitions and Q, defined in 
equation ((TJ, is the size of the Markov matrix. We use the identity 



{x + lf= f(x + l)'''' 




' 1 



ao,...,ai i=0 




E n -T L • (32) 



where J^i^'i = P- We remark that Q is the coefficient of in {x + 1)'^, 
whereas the coefficient of on the r.h.s. is precisely the number of ad- 
missible partitions defined in equation and is identical to the r.h.s of 
equation (p?T|l . 

In Table ^ the explicit example L = 15 and = 5 is worked out. We 
list the admissible partitions, calculate the corresponding order of degener- 
acy from formula ()28p and enumerate the corresponding multiplets by using 
equation P^jl . 

These results are invariant under 'particle - hole' exchange: the exclusion 
process with N particles jumping to the right can be mapped to a system 
with L — N particles jumping to the left after performing a particle - hole 
exchange. Of course the spectrum of the Markov matrix does not depend 
on the jumping direction. Thus, we know a priori that the Markov matrices 
of the TASEP with and L — N particles have the same spectrum. We 



Golinelli, Mallick — Spectral Degeneracies in TASEP 



12 



packages 


ao 


ai 


a2 


ag 


d{a) 


m{a) 


o o • • o 

• o o o o 
o • o o • 





5 








1 


243 


o o o • • 
o • o o • 


1 


3 


1 





1 


1620 


o o o o • 
• •ooo 
o • o o • 


2 


1 


2 





1 


810 


o • o o o 
• •ooo 
o • o o • 


2 


2 





1 


3 


90 


o • o o • 

o • o o o 
o • o o • 


3 





1 


1 


4 


15 



Table 1: An example of calculation of the degeneracies with L = 15 and 
= 5. Each row describes an admissible partition a = (oq, . . . ,ai) where 
tti counts the packages with i chosen roots. In the first column of this table, 
we have drawn one example of a choice set corresponding to a that selects 
N = 5 roots amongst the p = 5 packages of / = 3 roots each. The eigenvalues 
corresponding to a form m{a) multiplets of order of degeneracy d{a), where 
d{a) and m(a) are given by equations (j2Hll2ni)- The number of singlets {d = 1) 
is obtained by summing the first three contributions, so 2673 singlets. 

now verify this symmetry on the formulae derived above. Denoting with 
a 'tilde' the quantities for the model with N = L — N particles, we find 
that p = p (because gcd(L, A^) = gcd(L — N,N)), I = I and n = I — n. By a 
particle - hole exchange, a partition a = (ao, . . . , a^) is transformed to d where 
dj = ai^i. Hence, according to equations PH| I29 |1 . we obtain d{d) = d{a) and 
m{d) = m{a), i.e., degeneracies are indeed invariant by the particle - hole 
exchange. 

6 Discussion 

In this section, we draw some consequences of equations ()26II30|) . 

We first verify that these equations always predict the existence of singlets 
(i.e., isolated eigenvalues with degeneracy d = 1). In equation p8|l . we see 
that d = 1 ii and only if the partition has no empty package (oq = 0) or no 
full package {ai = 0) or if both ao = ai = 0. For instance, the stationary 
eigenvalue is always a singlet thanks to Perron- Frobenius theorem: in fact, 
the choice set for the stationary state is c = {1,...,A^} and hence each 
package Vk has n = N/p selected elements. The corresponding partition is 
thus given by = for i ^ n and a„ = p: this implies c? = 1, as expected. 

In order to obtain degenerate eigenvalues i.e., admissible partitions with 
d> 2, we must have Oq 7^ and a/ 7^ 0, i.e., at least one empty package and 
one full package. According to equation (j^ . the existence of degeneracies 
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L 




771,(1) 


m(2) 


777,(6) 


m(20) 


777(70) 


2 


1 


2 










4 


2 


4 


1 








6 


3 


8 


6 








8 


4 


16 


24 


1 






10 


5 


32 


80 


10 






12 


6 


64 


240 


60 


1 




14 


7 


128 


672 


280 


14 




16 


8 


256 


1792 


1120 


112 


1 


18 


9 


512 


4608 


4032 


672 


18 



Table 2: Spectral degeneracies in the TASEP at filling p = 1/2: L is the size 
of the lattice, N the number of particles; the other columns give m{d) the 
number of multiplets with degeneracy d. 



p 


L 




777(1) 


777(2) 


m(3) 


777(4) 


777(5) 


777(15) 


1/3 


9 

12 


3 
4 


81 
459 




1 
12 










15 


5 


2673 




90 


15 








18 


6 


15849 




540 


270 




1 




21 


7 


95175 




2835 


2835 


189 


21 


1/4 


16 
20 


4 
5 


1816 
15424 






1 
20 








24 


6 


133456 






240 


36 




1/5 


25 


5 


53125 








1 




2/5 


15 


6 


4975 


15 











Table 3: Examples of spectral degeneracies in the TASEP at filling p 7^ 1/2; 
L is the size of the lattice, N the number of particles; the other columns give 
m{d) the number of multiplets with degeneracy d. 

is given by the condition 1<N<L — I or equivalently by 

L<pN <pL-L (33) 

where p = gcd(L, A^). Some numerical examples are given in Tables [21 and El 
We now analyse the TASEP at a fixed value of the filling p = n/l, 77 and 
/ being mutually prime. The two integers (L, A^) are parameterised by the 
single number p defined in equation (j?T|) . According to the condition (jH^ . 
degeneracies appear when 
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This condition can always be fulfilled when < p < 1. Moreover, when the 
system size L and the number of particles increase with a given rational p, 
higher and higher orders of degeneracies appear and become generic in the 
large p (or L) limit. The admissible partition a = (ao, . . . ,ai) that maximises 
uj{a) = d{a) m{a) (i.e., the total number of eigenvalues associated with a) is 
given by 

a.~pQp^(l-py-V (35) 

The corresponding order of degeneracy increases exponentially with the size 
L as d (X where 

D=n+-\ fl + -j with u = p\ V = (I - p)K (36) 

For example, Z^i/a = 2^/'^ ^ 1.189 for p = 1/2; Dys = (9^/2^)^/^'^ ^ 1.040 
for p = 1/3, etc... D converges rapidly to 1 when the denominator / of p 
grows. That explains why, in numerical studies of systems of limited size, 
degeneracies are found only when I is rather small, i.e., when p is a 'simple' 
fraction. 

Similarly, we can also determine the admissible partition a that maximises 
m(a), the number of multiplets. This optimal partition has ai = when 
p < 1/2, and = when p > 1/2. In either case, this corresponds to 
d = 1. Thus the most numerous multiplets for p ^ 1/2 are singlets. This 
result does not contradict equations IHBj) in which the product d{a) m{a) 
is maximised. For the special case, p = 1/2, the partition that maximises 
m{a) has both ao and ai different from and the corresponding order of 
degeneracy increases as d oc 2^/6 (Golinelh and Mallick 2004). 

Furthermore, for a given number of particles, we can search the values 
of L where degeneracies appear. Because of the particle - hole symmetry 
(i.e., N L — N), we need to consider only the case < L/2. Then, 
the condition (jH!^ becomes 2A^ < L < pN. Because p < N, this implies 
that L < N"^: only a finite number of L values are possible. In the dilute 
limit (when L becomes large and remains fixed), the TASEP has thus no 
degeneracy. 

All these results have been derived on the basis of the 'one-to-one' hypoth- 
esis stated at the beginning of section El It is therefore crucial to compare 
our formulae with numerical results. We have numerically diagonalised the 
Markov matrix of the TASEP for certain values of the parameters {L,N). 
We use the translation symmetry to split the matrix of size fl into L ma- 
trices of size about Q/L. The spectrum is then computed by using Lapack 
library (Anderson et al 1999) and degeneracies are counted; details about this 
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procedure are given in (Golinelli and Mallick 2004). We have investigated 
systematically all the systems (L, A^) with L < 19. For L > 20, we have stud- 
ied the systems (L, N) in which non-trivial degeneracies are predicted and 
in which the size of the diagonalised matrix remains less than 6000. Results 
are given in Tabled for p = 1/2 and in Table El for other values of p. The 
numerical results are in perfect agreement with our analytical predictions, 
equations (|28IIH()|1 . 

The degeneracies for systems much larger than those listed in Tables |2] 
and El can be calculated from the formulae ()28IIH0j) for systems with sev- 
eral hundred sites and particles. However, the full numerical diagonali- 
sation of the Markov matrix consumes a computer time of the order of 

oc (p^(l — py^^)~^^. Such a fast growth limits the comparison between 
numerical diagonalisations and the exact formulae. 

7 Conclusion 

The spectrum of the Markov matrix of the TASEP on a one-dimensional 
periodic lattice has a rich structure with multiplets having degeneracies of 
high order. This structure depends on the filling fraction and presents arith- 
metical properties related to some particular partitions of the total number 
of particles. We have derived analytical formulae for the spectral degenera- 
cies by analysing the Bethe equations for the TASEP. These predictions have 
been verified by numerical calculations and we conjecture that the formulae 
we propose are exact. Our arguments are based on a 'one-to-one hypothesis' 
which is stronger than assuming the completeness of the Bethe Ansatz. Al- 
though the agreement between numerical results and analytical predictions is 
a strong argument in favour of the one-to-one hypothesis, it is possible that 
this hypothesis is not satisfied but that a weaker formulation, leading to the 
same spectral structure, holds good. We plan to study the completeness of 
the Bethe Ansatz and this one-to-one hypothesis more precisely in a future 
work. 

Our derivation of the spectral structure from the Bethe equations is rather 
indirect. In fact, the presence of such high orders of degeneracies is a strong 
evidence for hidden symmetries in the model. In other words, the TASEP 
should be invariant under a group operating on the configuration space such 
that each multiplet of a given order of degeneracy is an irreducible repre- 
sentation of this group. The orders of degeneracies would then classify the 
dimensions of irreducible representations and the number of multiplets of a 
given order would represent the multiplicity of a given irreducible represen- 
tation in the global representation over the configuration space. By using 
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the techniques of the algebraic Bethe Ansatz, we have constructed a series of 
nonlocal operators that act on multiple particles and that commute with the 
Markov matrix. We hope that a study of the representations of these oper- 
ators will allow us to understand the structure of the spectrum in a purely 
algebraic manner without having to analyse the solutions of the Bethe equa- 
tions. 

The TASEP evolution can be generahzed by introducing a fugacity pa- 
rameter A, i.e., by multiplying each non-diagonal term of the Markov matrix 

by a factor A. This parameter has been used to calculate the large deviation 
function of the total particle displacement (Derrida and Lcbowitz 1999). We 
verified that the spectral degeneracies for the 'TASEP + fugacity' model 
are the same as those found for TASEP. This is not surprising because the 
properties of the Bethe equations are not altered by introducing a fugacity 
and the arguments given in sections 4 and 5 can be generalized without any 
difficulty. However, in the case of the partially asymmetric exclusion process, 
in which particles can jump in both directions, the spectrum has a much 
simpler structure: eigenvalues are either singlets or doublets. Lastly, for the 
TASEP on an open lattice, the spectrum is made only of singlets. 
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Appendix: Layout of the solutions of the Bethe 
equations 

In this Appendix, we explain how to label the L solutions {Zi, . . . , Zl) of the 
polynomial equation of degree L 

(i-z)^(i + z)^-^ = y, (37) 

(0 < A'" < L), in such a way that each Zk{Y) is an analytic function of the 
parameter Y in the complex plane with a branch cut along the real semi-axis 
[0,+oo). 

A non-zero complex number Y can be written in a unique way as 

y = e*'^ with < < 2% , (38) 

r being a positive real number. This determination of the argument has a 
branch cut along [0, -|-oo). For a given value of r, the complex numbers Zk 
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Figure 1: Labelling the roots of the Bethe equations. Here L = 15, = 6, 
(f) = 7r/2 and r/rc = 0.8,1,1.2 (see text). The (resp. blue, black and red) 
continuous curves are the corresponding Cassini ovals. When r is fixed and 
(j) varies from to 2ti, each slips counterclockwise along a part of the 
Cassini oval. When cj) is fixed and r varies from oo to 0, each travels 
along a dashed curve from oo to points +1 or -1. 

belong to the generalized Cassini oval defined by 

\Z -l\f\Z + = r (39) 

where p = N/ L is the filling of the system. As shown in Fig. Q the topology 
of the Cassini oval depends on the value of r with a critical value 

r, = 2pP{l-py-^: (40) 

• for r < Tc, the curve consists of two disjoint ovals with N solutions on 
the oval surrounding +1 and L — N solutions on the oval surrounding 
-1. 

• for r = Tc, the curve is a deformed Bernoulli lemniscate with a double 
point at Zc = 1 — 2p. 

• for r > Tc, the curve is a single loop with L solutions. 

The Cassini ovals are symmetrical only if p = 1/2. 

In order to label the solutions, we start by considering the limit r oo 
for a given 0. Equation ()37|) then becomes 

Z^ ~ exp[i{(f) - Nn)] . (41) 
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The solutions are labelled by 



r exp 



- A^TT + 2(A; - l)7r] 



with = 1, 



(42) 



In other words, the Zk are regularly distributed along a large circle of radius 
r with 

^ < argZi < ... < argZi < + 2tt. (43) 

This labelling, obtained for large r, is extended by analytic continuation 
to all values of r, keeping (p fixed. The loci of the Z^ are drawn in Fig. ^ 
(dashed curves): they are orthogonal to the Cassini ovals. A singularity 
appears along the branch = because Zi and Z^^i collapse into each 
other at the double point Zc when r = rc\ we circumvent it by choosing 
= 0+. 

With this labelling, the solutions are ordered along the Cassini ovals. 
Moreover, when r < r^, the solutions {Zi, . . . , Zjq) group together on the 
right oval and {Z^+i, • • • , Zl) on the left oval. 
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